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ABSTRACT 



Context. In a recent work, we tried to obtain a calibration of the two components of the pre-main sequence binary system RS Cha by means of 
theoretical stellar models. We found that the only way to reproduce the observational parameters of RS Cha with standard stellar models is to 
decrease the initial abundances of carbon and nitrogen derived from the GN93 solar mixture of heavy elements by a few tenths of dex . 
Aims. In this work, we aim to reproduce the observational properties of the RS Cha stars with stellar evolution models based on the new 
AGS05 solar mixture recently derived from a three-dimensional solar model atmosphere. The AGS05 mixture is depleted in carbon, nitrogen 
and oxygen with respect to the GN93 mixture. 

Methods. We calculated new stellar models of the RS Cha components using the AGS05 mixture and appropriate opacity tables. We sought 
models that simultaneously satisfy the observations of the two components (masses, radii, luminosities, effective temperatures and metallicity). 
Results. We find that it is possible to reproduce the observational data of the RS Cha stars with AGS05 models based on standard input physics. 
From these models, the initial helium content of the system is K ~ 0.255 and its age is ~ 9.13 ± 0.12 Myr 
Conclusions. 
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1. Introduction 

Valuable tests of stellar physics are provided by the study of 
stars that are members of binary systems. When the global pa- 
rameters of both components of such systems are determined 
accurately from observations, the calibration of the system by 
means of stellar evolution models may yield the fundamental 
properties of its members, like the age or initial helium con- 
tent, as well as constraints on the input physics of the models. 

We recently studied the double-lined eclipsing binary sys- 
tem RS Cha, the components of which are both on the pre-main 
sequence (PMS). We obtained spectroscopic data that allowed 
us to determine accurate values of the masses and radii of the 



comp onents and the metallicity of the system (lAlecian et al 



2005h . We then complemented these parameters with previous 
determinations of the components' effective temperatures and 
luminosities (Table ^ and we attempted to calibrate the sys- 
tem, i.e. we tried to fit the observational data by means of stel- 
lar evolution models that we calculated for the two stars for 
a common value of the age and initial chemical composition 
( lAlecian et al. I l2006l lereafter paper I). We showed that stan- 



dard stellar models based on the iGrevesse & Noelsl (Il993h so- 
lar mixture (hereafter GN93) are unable to satisfy the observa- 
tional constraints for both stars even with reasonable changes in 
the model inputs. On the other hand, we found that both com- 
ponents can be fitted by standard stellar models of the same 
age and initial composition provided we decrease the C and N 
abundances of the GN93 solar mixture by a few tenths of dex. 

Recently, t he sol a r abun dances have been revised by 



Asplund et al] ( |2004 l2005allbl) on the basis of a time 



dependent, three-dimensional hydrody namical model of the 
solar atmosphere. In their analysis, lAsplund et al.l (12004 
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2005allbl) have considered departures from LTE and have used 
improved atomic and molecular data. T he new solar mixture 
which is listed in lAsplund et al. I (l2005bh (hereafter AGS05) is 
depleted in C, N, O with respect to the GN93 mixture. This in 
turn leads to a decrease in the solar (Z/X) ratio (ratio of the 
mass of the heavy elements with respect to hydrogen). 

In this paper, we present new stellar models based on the 
AGS05 solar mixture and appropriate opacity tables and we 
discuss their ability to reproduce the observational constraints 
of the two components of the RS Cha PMS binary system. 
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Table 1. Fundamental parameters of RS Cha. P stands for pri - 
m ary and S for second ar y References: l:lAlecian et alj j2005b 
2: lRibaset al](l2000l) . 3: IClausen & Nordstroml ( ll980l) 





Primary 


Secondary 


References 


M/Mo 


1.89 ±0.01 


1.87 ±0.01 


1 


R/Ro 


2.15 ±0.06 


2.36 ± 0.06 


1 


reff(K) 


7638 ± 76 


7228 ± 72 


2 


logL/Lo 


1.15 ±0.09 


1.13 ±0.09 


L = AnR^aT*^ 


logg 


4.05 ± 0.06 


3.96 ± 0.06 


g = GM/R^ 


Piday) 


1.67 


1 


i(degrees) 


83.4 ± 0.3 


3 


[Fe/H] 


0.17 


±0.01 


1 



2. Physical ingredients and initial parameters of 
the stellar models 



We used the CESAM stellar evolution code (lMorellll997l) to 
compute pre-main sequence stellar evolutionary tracks of con- 
stant mass in the range 1.86 to 1.90 Mq coiTesponding to the 
range of masses of the RS Cha components (Table [U. We do 
not take into account the effects of rotation nor those of mag- 
netic fields, and we neglect microscopic diffusion of the chem- 
icals (see paper I). Each evolution is initialized with a homo- 
geneous, fully convective model in quasi-static contraction and 
we define the age of a star as the time elapsed since initializa- 
tion. 

We calculated two series of models, each based on a dif- 
ferent chemical mixture of heavy elements. Hereafter we refer 
to GN93 models, the models based on the Grevesse & Noels 



( 119931) solar mixture and to AGS05 mode ls, the models based 
on the revised solar mixture obtained by .Asplund et al.l ( 12004 , 
2005allbl) . The GN93 and AGS05 solar mixtures differ. The 
abundances of C, N, O of the AGS05 mixture are lower which 
in turn leads to a large decrease of the solar (Z/X) ratio: 
(Z/X)q = 0.0245 for the GN93 mixture and (Z/X)q = 0.0171 
for the AGS05 mixture. 

Convection is treated using the clas sical mixing-length the- 
ory of convection ( Bohm- Vitensd 1958 ). The mixing-length pa- 
rameter value is taken to be equal to 1 .62Hp (Hp is the pres- 
sure scale-height). This value was obtained for the calibration 
of a solar model in luminosity and radius with the same input 
physics and the GN93 abundances. However the calibration of 
the solar model with the AGS05 mixture would change the so- 
lar mixing-length by less than Q.lHp. Furthermore as both stars 
of RS Cha enclose a very small convective core and a radiative 
envelope, the change in the value of the mixing length parame- 
ter does not affect their position in the HR diagram. 

We used the OPAL equation of state jRogers et al. 1996 ) 
and opacities calculated for the mixture used in the mod- 
els. In the interior we used the OPAL opacities tables 



dlglesias & Rogerslll996h for either the GN93 or the AGS05 
mixture (we generated these latter on the OPAL Web site) 
while a t low temperatures (i.e. for T < 10^ K), we used ei- 
ther the lAlexander & F erguson (1994j) o pacitie s -provided for 



The species entering the nuclear reaction network are : 'H, 
^He, ^He, '^c, i^C '^n^ i5n^ isq, '^O. We considered ^H, ^Li 
and 'Be to be in equilibrium. The r iuclear reaction rat es are 
taken from the NACRE compilation dAngulo et al.lll999h . 

In our models, we take Z/X = 10^^''^^\Z/X)q, where 
[Fe/H] is the observed value given in Table [T] and {Z/X)q is 
either the value of GN93 or the value of AGS05. We then de- 
rive the individual abundances of the heavy elements from the 
chosen solar mixture (GN9 3 or AGS05). 



We are aware that in lAlecian et al.l (120051) we obtained 
the metallicity of th e system from a one dimensional ATLAS 
model atmosphere ( lKuruczlll993 ) calculated with the GN93 
old solar mixture. To be consistent, to calculate internal struc- 
ture models based on the new solar mixture, we should 
have redetermined the [Fe/H] -value of RS Cha from a three- 
dimensional model atmosphere calculated with the AGS05 
mixture. However 3-D model atmospheres are still not avail- 
able for A-stars which have finer convection zones and larger 
convection cells. On the other hand, it would be incorrect to 
use 1-D Kurucz models with the AGS05 abundance to reder- 
mine the metallicity of RS Cha. 



Alecian et a l. (2005) have determined the metallicity of RS 



the GN93 mixture- or the .Ferguson et al.l(l2005h opacities avail 
able for the AGS05 mixture. 



Cha from measurements of the Fe and Ca abundances. Since 
the abundances of Fe and Ca are very close in the GN93 and 
AGS05 mixtures, it seems to be reasonable to assume that an 
analysis of the spectra with 3-D model atmospheres would not 
much change the [Fe/H] -value of the system, which is a differ- 
ential value. Under this assumption, we find that the revision of 
(Z/X)q from 0.0245 (GN93 mixture) to 0.0171 (AGS05 mix- 
ture) leads to a decrease of the (Z/X) ratio of RS Cha from 
0.0362 to 0.0253. 

However for stars more massive than the Sun, work in 
progress seems to show that 1-D and 3-D model atmospheres 
yield non negligeable differences in [Fe/H] (Thevenin, private 
communication). To investigate the possibility that the 3-D de- 
termination leads to a value of [Fe /H] very different from the 1 - 
D result of Table[T] we have calculated models with [Fe/H] val- 
ues spanning the range 0.07-0.27 dex, that is ±0.10 dex around 
our mean value of 0.17 dex. 

The GN93 models presented below have been calculated 
with an initial helium mass fraction = 0.267 derived from 
the calibration of the solar model in luminosity and radius. 
However we have shown in paper I that a change of the ini- 
tial helium content does not change the behaviour of the GN93 
models that cannot reproduce the observations. The calibra- 
tion of the solar model in luminosity and radius of the solar 
model with the AGS05 abundances leads to an initial solar he- 
lium content that is lower by ~ 0.01 than that obtained with 
the GN93 abundances. In the new AGS05 models calculated 
for the present study, we have ajusted the initial helium content 
Y to reproduce the observations, which yields Yq - 0.255 (see 
Sect.[3]below). 

3. Comparison with observations for GN93 and 
AGS05 stellar models 

The RS Cha system is composed of two similar A-type stars 
in a PMS stage of evolution close to the ZAMS at the onset of 
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Fig. 1. Comparison of the GN93 models with observations. 
Crosses are the observed error bars in effective temperature and 
luminosity (Table [T]| for the primary (P) component (left) and 
for the secondary (S) component (right). Dot-dashed lines are 
evolutionary tracks for the GN93 models of different masses: 
lowest to upper tracks are for masses of 1.86, 1.87, 1.88, 1.89 
and 1 .90 Mq. Boxes superimposed on the tracks with a full line 
correspond to models that satisfy the observed constraints on 
masses and radii. 



CNO burning. Figures [T] and |2] are H-R diagrams showing the 
comparison of the observations of the RS Cha components with 
stellar evolutionary tracks for the two mixtures considered. 

Observed effective temperatures and luminosities as well 
as masses and radii of the RS Cha components are compared to 
the corresponding values obtained from stellar model, as fol- 
lows. For a given mass Mobs given in Table [T] a stellar model 
is evolved until its radius matches the observed radius /?obs- 
This provides the calculated luminosity and effective temper- 
ature (L, TeftOcaic of the considered star. The same is done for 
the extreme values of the observed masses and radii and this 
gives rise to error boxes in effective temperature and luminosity 
which are superimposed onto the tracks in Figs. [T] and |2] with 
full lines. In the same graph, crosses represent the error bars 
in luminosity and effective temperature for the primary (left) 
and secondary (right) components derived from observations 

(L, reif)obs- 

When both error boxes are located on their respective 
crosses, the corresponding stellar models reproduce the obser- 
vations. As it can be seen in Fig. [1] this does not occur when 
GN93 models are used. On the other hand. Fig. |2] shows that 
the AGS05 models calculated with Y - 0.255 well match the 
observations. We point out that the AGS05 models can repro- 
duce the observations without making any other change in the 
(standard) input physics of the models. 

As shown in paper I the GN93 models cannot reproduce 
the observations, even if we consider reasonable uncertainties 
in the physical inputs or if we change their initial helium con- 
tent and their metallicity within the error bars. The main prob- 
lem is that GN93 models predict that the luminosity of the pri- 
mary (Lp) is smaller than that of the secondary (Ls) which is 
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Fig. 2. Same as Fig.[T]for the AGS05 stellar models. The dots, 
diamonds and squares indicate the ages of the models. 



not supported by the observations. We also showed in paper I 
that, in order to invert the luminosity ratio Lp/Ls, we need to 
delay the luminosity drop of the tracks toward the main se- 
quence, induced by the onset of a convective core in the stars. 
This can only be obtained by decreasing the abundances of C 
and N. Indeed, as we remove carbon and nitrogen, the nuclear 
reactions involved in the CNO cycle become less efficient, and 
therefore the nuclear energy (e) produced in the center of the 
star decreases. We showed that the decrease of e leads to a 
dominant decrease of the ratio L/m, L and m being the local 
luminosity and mass, compared to a negligible increase of the 
ratio p-, where k, p and T are the local opacity, pressure and 
temperature. The radiative gradient Vraj, proportional to fjfr. 
is therefore reduced, and the onset of the convective core is 
delayed. The AGS05 models therefore match the observations 
because the AGS05 mixture is depleted in C, N and O with 
respect to the GN93 mixture. 

In order to ensure that the agreement would be kept even if 
the [Fe/H] value of RS Cha was to be modified by an anal- 
ysis of the spectra based on three-dimensional model atmo- 
spheres, we calculated AGS05 models with [Fe/H]-values of 
0.07 and 0.27 dex. We find that these models are still able to 
reproduce the observations after an adjustment of their initial 
helium abundance. 

The calibration of the RS Cha system based on the AGS05 
stellar models and setting the same age for both components 
leads to an age of 9. 13 + 0. 12 Myr for this PMS system, which 
is similar to the age we obtained in paper I from models where 
the C and N ratios had been decreased in the GN93 mixture. 
This age is the intersection of the age ranges given by both er- 
ror boxes. The error bare is therefore very small, however it 
does not take into account the uncertainties of the physical in- 
put parameters such as the mixing length or the overshooting 
parameters. Recent studies show that changing these parame- 
ters in a reasonable range may lead to a modification of the age 
determination u p to +20 % depending on the star considered 
( Lebretonll2005 ). A reasonable error bar on our age determina- 
tion would therefore be +2 Myr. 
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4. Conclusions 

We close here a series of papers devoted to the observational 
and theoretical study o f the two component s of the PMS bi- 
nary system RS Cha. In I Alecian et alj (120051) . we revisited the 
fundamental parameters of the system and we determined its 
metallicity for t he first time from spe ctroscopic observations. In 
a second paper (lAlecian et al we modeled the two stars 

of the system and have shown that stellar models built with 
standard input physics and the standard GN93 sola r chemical 
mixture of heavy elements (l Orevesse & Noelsll993 ) cannot re- 
produce the observations. We found that the only way to repro- 
duce the observations is to decrease the C and N abundances by 
0.6 and 0.5 dex respectively with respect to the GN93 mixture. 

Here, we have obtained a satisfactory calibration of the RS 
Cha PMS binary system with stellar evolution models based 
on standard input physics and on the new AGS05 solar mixture 
recently derived by Asplund et al. (2004, 2005a b). This is due 
to the fact that carbon, nitrogen and oxygen are depleted in the 
AGS05 solar mixture with respect to their values in the GN93 
mixture. 

To reproduce the observations, there are two alternatives: 

- Th e solar mixture of t he he avy elements corresponds to 
the Grevesse & Noelsl (Il993 h solar mixture. In this case, 
the abundance ratios of the chemical species in RS Cha 
cannot be solar an d are typical of stars in young clusters 
JPaflon et al.ll2004 . 

- The sol ar mixture of the heavy elem ents corresponds to the 
revised lAsplund et al.l(l2004l2005allbl) solar mixture. In this 
case the abundance ratios of the species in RS Cha are com- 
patible with a solar mixture. 



As pointed out by lAsplund et al. I (l2005bh the new AGS05 
mixture corresponds to C and O abundances which are in better 
agreement with those measured in the local interstellar medium 
and nearby B stars. However they also make the agreement 
between the solar mo del and hehoseismic observations much 
worse (see for instance ^Basu & Antiall2004 ). It would therefore 
be interesting to bring independent arguments in favor of (or 
against) the AGS05 mixture. 

The available observations of RS Cha and of other mem- 
bers of the T] Cha young stellar cluster cannot favor one of the 
alterna tives given above. Recent studies of the 77 Cha cluster 
I et all J: 



by lLyo 1 



(12004 1200 6) identify 18 stars as members of the 



cluster, among which 15 are K-M type stars. A possibility to 
discriminate models calculated with the two mixtures would 
be to compare the corresponding isoc hrones with obser vations 
in the color-magnitude diagrams of iLvo et al. I (l2004 I2OO6,) . 
We checked that the isochrones calculated with CESAM for 
the GN93 m ixture and a solar metallicity are close to those of 
Siess et al.l (12000) . Furthermore we have estimated the age of 
the 77 Cha cluster using RS Cha and the 5 m ore massive stars of 
the low-mass star sample of lLvo et all (12004 ) (RECX 1, 4, 7, 10 
and 11), for both mixtures : GN93 and AGS05. In both cases, 
we found a large ag e range co nsistent with the one found by 
Lawson et all (l200lh . from the ID' Antona & Mazzit elli ( 1997, 
tracks, and by ILvo et al.l (|2004 from the 



Siess et al., 



1998 

(I2OOOI) tracks. However for RECX 1, 7 and 10, being binaries 



(ILvo et al.ll2004 . we can reduce the age range by taking into 
account solely the single stars RECX 4 and 1 1 . The age range 
obtained is consistent with our age determination of RS Cha 
whatever the mixture considered. Therefore we are not able to 
test the different solar mixtures with the low-mass stars of the 
77 Cha cluster with the isochrones fitting procedure. 

Another possibility would be to study the chemical com- 
position of the 77 Cha cluster by determining the carbon, nitro- 
gen or oxygen abundances in the late-type stars of the cluster 
and to compare them to both solar mixtures. These abundances 
have not been determined yet, and spectroscopic observations 
of these stars would be helpful to support one of these alterna- 
tives. 

In both cases, the models predict that the initial helium 
abundance of the system is subsolar and that its age is of about 
9 ±2 Myr, w hich is consistent with the recent age determination 
of lLvo et al.i(l2004 120061) . 

Acknowledgements. We are very grateful to Frederic Thevenin for 
fruitful discussions. 
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